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Abstract 

We study structural properties of the g-color Potts field theory which, for real values of q, 
describes the scaling limit of the random cluster model. We show that the number of in- 
dependent n-point Potts spin correlators coincides with that of independent n-point cluster 
connectivities and is given by generalized Bell numbers. Only a subset of these spin correla- 
tors enters the determination of the Potts magnetic properties for q integer. The structure 
of the operator product expansion of the spin fields for generic q is also identified. For the 
two-dimensional case, we analyze the duality relation between spin and kink field correlators, 
both for the bulk and boundary cases, obtaining in particular a sum rule for the kink-kink 
elastic scattering amplitudes. 



1 Introduction 



There are cases in which the notion of spontaneous symmetry breaking which is commonly used 
to describe ferromagnetic phase transitions admits an extension which allows the description of 
phase transitions of a different nature. The case of interest in this paper is that of the g-state 
Potts model [H [2] , whose Hamiltonian 

= -J^^s{x),s{y) , S{x) = l,...,q, (1) 

{x,y) 

is invariant under the group Sq of permutations of the values of the variables s{x) located at 
each site of a lattice C. The q values of s(x) can conveniently be thought as q different colors 
and, clearly, the Potts ferromagnet is well defined for an integer number of colors. It has been 
known for a long time, however, that the partition function X^j^^^,)! can be rewritten, up to 
an inessential constant, as [3] 

Z=^p->'il-pr>'q^% (2) 

where ^ is a graph obtained putting rib bonds on the lattice C, each one with weight p = \—e~^ € 
[0, 1] (rife is the number of absent bonds in C). The connected components of including isolated 
sites, are called Fortuin-Kasteleyn (FK) clusters. In the formulation ([2]), which defines the so 
called random cluster model, the probability measure for the graph Q depends on q through the 
factor being the number of clusters in ^, and is well defined for any real positive q. In 

the thermodynamic limit the random cluster model undergoes, for q continuous, a percolative 
phase transition associated to the appearance, for p larger than a (^-dependent critical value 
of a non-zero probability of finding an infinite cluster. The transition is first order for q larger 
than a dimensionality-dependent value qc^ and second order otherwise; in particular, the limit 
g — )• 1, which eliminates the factor q^" in ([2D, describes ordinary percolation. When q is an 
integer larger than 1, the percolative transition of the random cluster model contains, in a sense 
to be clarified below, the ferromagnetic transition of the Potts model. 

For q < qc, when a scaling limit exists, the problem admits a field theoretical formulation. 
There must be a field theory, which we call Potts field theory, that describes the scaling limit of 
the random cluster model for q real, as well as that of the Potts ferromagnet for q integer. This 
theory has the Potts spin fields as fundamental fields (the FK mapping relates Potts spin corre- 
lators and connectivities for the FK clusters) and is characterized by S'g-invariance. The obvious 
question of the meaning of Sq symmetry for q non-integer arose at least since the e expansion 
treatment of [1], and appears to admit a general answer: although one starts from expressions 
which are formally defined only for q integer, formal use of the symmetry unambiguously leads 
to final expressions containing g as a parameter which can be taken continuous. 

The two-dimensional case allows for the most advanced, non-perturbative results. The Potts 
field theory is integrable and the underlying scattering theory was exactly solved in [5] for 
continuous q < qc = ^- The scattering solution was then used in [6l [7] to determine two-point 
correlators and universal combinations of critical amplitudes for the Potts model and percolation. 
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In spite of these results, even in two dimensions the nature of Potts field theory for real values of 
q needs to be better understood. Indeed, already in the critical case, although exact formulae for 
the scaling dimensions and the central charge as functions of q have been known for a long time 
[HI [9] , and results for boundary correlators have been obtained since the work of [10] , a complete 
conformal field theory description is only available for the magnetic properties at g = 2, 3, 4. The 
essential problem is readily stated: since the Potts spin field has q — 1 independent components, 
no ordinary characterization of the space of fields appears possible for q non-integer. Recently 
[llj . we showed how the problem of non- integer field multiplicities can be dealt with for three- 
point functions, and arrived in this way to exact predictions for operator product expansion 
(OPE) coefficients which have already been confirmed by high-precision Monte Carlo simulations 
for the case of ordinary percolation [12\ [13] . 

In this paper we investigate, for q continuous, some structural properties of the Potts field 
theory as a theory characterized by Sq invariance under color permutations and able to describe 
the scaling limit of the random cluster model. We first of all observe that the issue of the content 
of the theory is better addressed, in any dimension, focusing on linearly independent correlation 
functions rather than on field multiplicities. For this purpose one needs to have in mind the 
relation of spin correlators with cluster connectivities for q real, rather than with magnetic 
properties for q integer. Just to make an example already considered in [11) . the correlator of 
three spin fields with the same color is proportional to the probability that three points are in 
the same FK cluster. This probability is well defined and non-vanishing for continuous values 
of q in the random cluster model, in particular for the case g' = 2 in which the spin correlator 
has a zero enforcing Ising spin reversal symmetry; stripped of a trivial factor q — 2, this spin 
correlator enters the description of cluster connectivity at q = 2 as for generic real values of 
q. Similarly, a number nc < n of different colors enters the generic n-point spin correlator: a 
correlator with ric > q has no role in the description of the Potts ferromagnet, but enters the 
determination of cluster connectivities in the random cluster model. One then realizes that the 
dimensionality F„ of the space of linearly independent n-point spin correlators for real values 
of q is actually (/-independent and must coincide with that of the space of linearly independent 
n-point connectivities. We show that this is indeed the case and that Fn coincides with the 
numbei0 of partitions of a set of n elements into subsets each containing more than one element; 
the relation between spin correlators and cluster connectivities is also given and written down 
explicitly up to n = 4. Only a number Mn{q), smaller than Fn for n large enough, of independent 
spin correlators enters the determination of the magnetic properties at o integer, making clear 
that the magnetic theory is embedded into the larger percolative theorjo. 

Once the relevant correlation functions have been identified, an essential tool for their study 
is the OPE of the spin fields. Again, the existence of such an object for real values of q is made 
a priori not obvious by the badly defined multiplicity of the fields. We show, however, that its 
structure can be very naturally identified (equations (US]) and (pU]) below). 

^We consider the symmetric phase, i.e. the case p < pc. 

^See |14l 115) for detailed studies of this fact in the Ising model. 
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The additional property we study in this paper, duahty, is specific of the two-dimensional 
case. It is well known [Tl |2] that spin correlators computed in the symmetric phase of the square 
lattice Potts model coincide with disorder correlators computed in the spontaneously broken 
phase. Here we study duality directly in the continuum, for real values of q, and with the main 
purpose of clarifying the role of kink fields. These are the fields that in any two-dimensional 
field theory with a discrete internal symmetry create the kink excitations interpolating between 
two degenerate vacua of the spontaneously broken phase; in general, they are linearly related to 
the usual disorder fields, which are mirror images of the spin fields, in the sense that they carry 
the same representation of the symmetry. In the Potts field theory, however, it appears that in 
many respects the use of kink fields provides a simpler way of dealing with the symmetry for real 
values of q. The duality between spin and kink field n-point correlators is a non-trivial problem 
that we study in detail up to n = 4, both for bulk and boundary correlations. The problem 
of correlation functions for points located on the boundary of a simply connected domain is 
simplified by topological constraints which, in the limit in which the boundary is moved to 
infinity, also account for non-trivial relations among kink scattering amplitudes. 

The paper is organized as follows. In the next section we investigate cluster connectivities 
and Potts spin correlators, their multiplicity and the relation between them. In section 3 we 
analyze OPE's and obtain in particular that for the Potts spin fields for real values of q. Duality 
between spin and kink field correlators in two dimensions is studied in general in section 4 and 
specialized to boundary correlations in section 5. Few final remarks are given in section 6 and 
three appendices complete the paper. 

2 Counting correlation functions 
2.1 Cluster connectivities 

Correlations within the random cluster model ([2]) are expressed by the connectivity functions 
giving the probability that n points xi,...,Xn fall into a given FK cluster configuration. In 
order to define the connectivities we associate to a point Xi a label Oj, with the convention that 
two points Xi and Xj belong to the same cluster if Oj = aj, and to different clusters otherwise. 
We then use the notation -Pai...a„(3^i) •■■^Xn) for the generic n-point connectivity function, within 
the phase in which there is no infinite cluster, i.e. for p < Pc- The total number of functions 
Pai...a„{xi, ■ ■ ■ ,Xn) is the number i?„ of possible partitions (clusterizations) of the n pointjfl; 
these Bn functions sum to one and form a set that we call C^"'\ 

It is not difficult to realize that the elements of C'-"-' can be rewritten as linear combinations 
of "basic" A;-point connectivities, with k = 2, . . . , n; we call the number of basic /c-point 
connectivities, and V^'^^ C C^^^ the set they form. There is a simple procedure to build P^"-) 
given the ■pC'^'s for /c = 2, . . . , n - 1. Let us start with n = 2: C^^) = {Paa,Pab}, but Paa + 
Pab = 1 implies F2 = 1 and we can choose "P*-^-* = {Paa}- Consider now n = 3: C(3) = 

^The B„'s are known as Bell numbers and are discussed in Appendix A. 
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aaa 


+ Paab 


= PaaO 


= Paa{xi,X2 


P 

^ aaa 


+ Paba 


= PaOa 


= Paa {■^l 1 -^3 


P 

^ aaa 


~l~ Pbaa 


= Poaa 


= Paaix2,X3 



{Paaa, Paab, Paba, Pbaa, Pabc}, however starting from Paaa and Summing over the inequivalent 
configurations of the last point we have Paaa + Paab = PaaO, where PaaO = Paaixi,X2) belongs, 
we chose it on purpose, to V^'^K We observe that only one linear combination of connectivities 
in C^^^ reproduces, taking into account the coordinate dependence, the two-point connectivity 
in p(^) ; we obtain then the sum rules 

(3) 
(4) 
(5) 

Paaa ~l~ Paab ~l~ Paba ~l~ Pbaa ~l~ Pabc — !• (6) 

This system of equations exhausts all the possible linear relations among the elements of C^^^; 
a reduction to a non-basic two-point connectivity (i.e. not belonging to 'P^'^\ for example Pabo) 
will indeed produce an equation which is a linear combination of those above. It follows, in 
particular, that the five elements of C^^^ can be written in terms of Paaa, PaaO, Poaa, PaOa, so that 
F3 = 1 and we can choose V^^^ = {Paaa}- 

In general, suppose we have chosen the Fk basic connectivities in V^'^^ for k = 2, . . . ,n — 1. 
Given -Pai...a„ S C^"'\ we can fix k of its n indices according to an element of V^^^; summing over 
the inequivalent configurations of the remaining n — k indices we obtain a linear relation for the 
connectivities of C^'"\ We can do this for each of the Fj. elements in V^'^^ and for choices 
of k indices among n indices. The number of independent sum rules for the elements of C*-""^ is 
then 

with the last term accounting for the fact that the n-point connectivities sum to one. The 
number Fn of elements of "P^") is the minimum number of connectivities in C*-"^ needed to solve 
the linear system of £"„ equations in i?„ unknowns, i.e. 

Defining Fq = 1 and knowing that Fi = 0, we rewrite ([8]) as 

Pn = J2(t)^^^ Vn>0. (9) 

fc=0 ^ ^ 

We show in Appendix A that ([9]) implies that Fn is the number of partitions of a set of n 
points into subsets containing at least two points. We list in Table [1] the first few Bn and Fn- 
The combinatorial interpretation of the FnS suggests that a natural choice for the set of 
linearly independent n-point connectivities is to consider clusterizations with no isolated points, 

i.e. = {p^^}^ ^(3) ^ p(4) = {Paaaa,Paabb,Pabab,Pabba}, and SO On. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Bn 


1 


2 


5 


15 


52 


203 


877 


4140 


21147 


115975 


Fn 





1 


1 


4 


11 


41 


162 


715 


3425 


17722 


Mn{2) 





1 





1 





1 





1 





1 


M„(3) 





1 


1 


3 


5 


11 


21 


43 


85 


171 







1 


1 


4 


10 


31 


91 


274 


820 


2461 



Table 1: The Bell numbers Bn give the number of partitions of n points. The number Fn 
of hnearly independent n-point spin correlators (|12|) coincides with the number of partitions 
on n points into subsets containing at least two points. A number Mn{q) of these correlators 
determines the ?z-point magnetic correlations in the g-state Potts ferromagnet. 



2.2 Spin correlators 

As we will see in a moment, it follows from the FK mapping that the Potts spin correlators can 
be expressed as linear combinations of the cluster connectivities. Consistency of this statement 
requires that the number of linearly independent spin correlators coincides with the number of 
linearly independent cluster connectivities. The spin variables of the Potts model are defined 
afl 

aa{x) = q6s(x),a - , a = l,...,q, (10) 

where s{x) is the color variable appearing in ([1]), and satisfy 

q 

^aUx)=0. (11) 

The expectation value (ctq) is the order parameter of the Potts transition, since it differs from 
zero only in the spontaneously broken phase. More generally we denote by 

Gai...an , . . . , Xn) = {'-^ai {Xi) ...aa,XXn))j<J, (12) 

the n-point spin correlators in the symmetric phase. We now show that, for q real parameter, 
the number of linearly independent functions (jl2p coincides with F„. 

Let the string (ai . . . a^) identify the correlator ([12]) . and suppose that is isolated within 
this string, i.e. it is not fixed to coincide with any other index within the string. We can then 
use pip to sum over at and obtain, exploiting permutational symmetry, a linear relation among 
correlators involving a string similar to the original one together with strings without isolated 
indices. The simplest example, 

= ^ Gal3 = Gaa + (? " l)^^^ , 7 / « ; (13) 

is sufficient to understand that (jlip produces meaningful equations also if q is non-integer, 
the only consequence being that some multiplicity factors in front of the correlators become 



*Our present normalization of aa{x) difTers by a factor q from tliat used in [TT] . 
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non-integer; also, the requirement 7 7^ a should imply q > 2, but in the sense of analytic 
continuation to real values of q (|13p is equally valid for q < 2. Similarly, starting with a string 
with m isolated indices a^^ , • • • , ; summing over a^^ and using permutational symmetry we 
generate a linear relation involving the original string with the m isolated indices , • • • , cx-k^ 
together with strings in which is no more isolated, i.e. with at mos1§| m — 1 isolated indices. 
Iterating this procedure, any correlation function (I12p with isolated indices can be written as 
linear combination of correlation functions without isolated indices. Recalling the meaning of 
Fn in terms of set partitions, we then see that there are at most Fn Sg-inequivalent, linearly 
independent correlation functions (jl2p . i.e. those without isolated indices. On the other hand, 
the constraint (llip does not generate any new linear relation if we start with a string which 
does not contain any isolated index. The number of linearly independent functions (jl2p is then 
exactly F„. 

Let us now detail the linear relation between cluster connectivities and spin correlators. The 
latter admit the FK clusters expansion, 

s{xi),...,s{xn) GQC i=l 

where the prime on the first sum means that sites belonging to the same FK cluster are forced 
to have the same color s. Notice that, if one of the points Xi is isolated from the others in a 
cluster of a given graph Q, then the sum over its colors gives zero due to ([TT]) : hence, consistently 
with our previous analysis, (jl4p receives a contribution only from partitions of the sites rr^'s into 
clusters containing at least two of these sites. If i-*ai...a„ (a^i, • • • ,Xn) is the probability of such a 
partition, then the number of distinct clusters will be m < n, and to any pair (xj,aj) we can 
associate one of the distinct letters ci, . . . , Cm chosen among the Oj's. The coefficient of Pai...a„ 
in the expansion (jl4p is ther^l 

^ E n i^^s^a. - 1) • • • E n i'i^^r.,.. - 1) ; (15) 

the notation Xi C c means that to the point Xi is associated the letter c (x belongs to the cluster 
c). For n = 2, 3 the dimensionality of correlation spaces is F2 = F3 = 1 and gives 

Gaa = qiPaa, (16) 

Gaaa = 9l92 Paaa , (17) 

where we introduced the notation 

qk = q-k. (18) 

The first relation with a matrix form appears at the four-point level {F4 = 4), for which (jl5p 



""We can have strings with m — 1 or m — 2 isolated indices. 

^The prefactor l/g™ ensures the correct probabihty measure for the graph Q. 
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leads to 



Gaaaa = Qliq - 3g + 3)Paaaa + Ql {Paabb + Pahba + Pabab), (19) 

GaaPP = (2(7 — 'i)Paaaa + Q^Paabb + Pabba + Pabab, (20) 

GaPPa = (2(7 — ^)Paaaa + Paabb + QlPabba + -fafeab; (21) 

GaPaP = (2(? — 3)Paaaa + -Paafeft + Pabba + Q^Pabab- (22) 



The last set of equations, as well as those one obtains for n > 4, can be inverted to express 
the connectivities in terms of the spin correlators, making clear that all the F„ inequivalent 
and independent spin correlators are necessary to determine the connectivities of the random 
cluster model. On the other hand, a number Mn{q) < Fn of these spin correlators determine 
the magnetic correlations in the Potts model at q integer. This is due to the fact that the 
spin correlators are themselves the magnetic observables, and that for q integer some of them 
vanish (e.g. Gaa...a{xi, ■ ■ ■ ,x„) at q = 2, n odd, see (fT7|) ). those involving more than q colors 
are meaningless in the magnetic context, and additional linear relations may hold at specific 
valued of q. The numerical sequences Mn{q) are determined in Appendix B for g = 2,3,4 (the 
case (7 = 2 is of course trivial) ; the first few values are given in Table [1] and plots are shown in 
Fig. [TJ It is interesting, in particular, to compare the large n behavior of the dimensionalities of 
the magnetic and percolative correlation spaces. Defining 

M„(g)"^'e"^«("), (24) 

(jlOOP and (jlOSp give Sq{n) = \og{q — 1), a result which can be consistently interpreted as a kind 
of entropy for the Potts spin a a and is expected to hold for any integer q > 1. Fn exhibits 
instead the super-exponential growth [16] 

Fn e'""(") , (25) 

lost losr n / 1 \ 

s(n) =logn-loglogra- 1+ f ^ +0 . (26) 

log n \ log n J 

2.3 Scaling limit and correlators of kink fields 

For q < qc, i.e. when the phase transition is continuous, the Potts field theory describes the 
scaling limit J —t- Jc of the Potts model, with the spin variables cra{x) playing the role of funda- 
mental fields {x is now a point in Euclidean space). In particular, the q degenerate ferromagnetic 

'^For example the relation SGaaaa = 2{Gappa + Gapap + Gaa^p) holds specifically at g = 3 and the system of 
equations (19)-(22) is no longer invertible. More generally, we expect that the Fn x F„ matrix T„{q) giving the 
spin correlators in terms of the "basic" connectivities (given explicitly by (|16|) . (|17|l and (|19|) - (|22|l for n = 2, 3,4) 
has determinant 

n-l 

det(T„) = q"" (g _ 1) J] (q - fc)'^'^(^-) , (23) 

fc=2 

with dn{k) = X]j=i S{n,j) — A4„{k), S{n,j) being the generalized Stirling numbers discussed in Appendix A, and 
a„ determined by the requirement that the total degree of the polynomial (|23|) is Dn = X]fc=i('^ ~ k)S{n, k), as 
follows examining (fTS)). It follows from (|93)l and (|106|l that d„(A;) = for A; > n. We checked ((23)l up to ?^ = 5. 
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Figure 1: Plots of the first 20 values of the sequences Fn (number of independent n-point spin 
correlators for q real) and Mn{q) (number of independent n-point spin correlators in the magnetic 
sector) for q = 3,4. 



ground states which the Potts model possesses above Jc correspond in the scaling limit to de- 
generate vacua of the Potts field theory. In the two-dimensional case the kinks interpolating 
between a vacuum with color a and one with color /3 are topologically stable and provide the 
elementary excitations of the spontaneously broken phas^; they are created by the kink fields 
IJ'a/six), which are non-local with respect to the spin fields aa- The products of kink fields are 
subject to the adjacency condition /ia/3(a;)///3^(y). Duality relates, in a way to be investigated 
in the next sections, the kink field correlators in the broken phase 

Gi3,...i3„{xi,...,Xn) = (/i/3i/32(xi)Ai/32/j3(2;2) . . .^/3„/3i(x„))j>j^ (27) 

to the spin correlators in the symmetric phase ()12p . Consistency of the duality relation requires 
that the number of S'g-inequivalent correlators ()27p coincides again with Fn- We now show that 
this is indeed the case. 

Consider to start with the string of kink fields M/3i/32/^/92/33 ■ ■ ■ I^PnPn+i ^^"^ associate to it n + 1 
points Pj,i = l,...,n + 1, ona line. Each point Pi has a color /3j which must differ from those of 
the adjacent points. Let us show first of all that the number C„+i of S'q-inequivalent colorations 
of the points Pi coincides with the Bell number B^- If the adjacency condition is relaxed, the 
number of inequivalent colorations of the n + 1 points is Bn+i- The string will consist of A; + 1 
substrings, each with a definite color different from those of the adjacent substrings, that we 
can think to separate by placing k domain walls between them; this can be done in (^) ways. 
The k + 1 substrings can be colored in C^+i 5q-inequivalent ways and we have 

5n+i = X^Qcfc+i Vn>0. (28) 
*See [17] for a lattice study of Potts kinks. 
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Since Ci = 1, the result Cn+i = Bn then follows from (j86|) by induction. 

The case we just discussed includes L^+i inequivalent colorations in which /3i = (the 
case ([27|) we are actually interested in) and On+i inequivalent colorations in which (3i ^ /^n.+i, 
i.e. Ln+i + On+i = Cn+i- On the other hand, if we start with n points having /3i 7^ /3„ 
and we add Pn+i with = Pi, the number of inequivalent colorations does not change, i.e. 
Ln+i = On- We then see that 

K+l + Ln+2 = Cn+1 =Bn, V?l > . (29) 

This relation, together with the initial condition Li = 1, can be used to generate the -L^'s from 
the BnS. Of course Ln+i is the number of inequivalent correlators (j27p we were looking for. 
Comparison of (j29p with (|9ip then leads to the final identification Ln+i = F^. 

3 Operator product expansions 

Generically the OPE of scalar fields Ai{x) with scaling dimension Xi takes the form 

^linr^ MxM,{x,) = Y: C^^^ , (30) 

m ^12 

where we include for simplicity only scalar fields in the r.h.s. and use the notation 
in the following we will replace (j30p by the symbolic notation 

Ai • Aj = ^ ^ Clj Am ■ (31) 

m 

The nature of the fields ^lapi^) naturally leads to the two-channel OPE 

Ma/3 • M/37 = ^a'^I + (1 ~ ^a'y){C^fla'y + • • •) ' (32) 

where the neutral channel = 7 contains the expansion 1 = 1 + C^s + • • • over Sg-invariant 
fields (identity I, energy e, and so on), and the charged channel a 7^ 7 the expansion over fi^^ 
and less relevant kink fields; and are simplified notations for the OPE coefficients, for 
which exact expressions for continuous q have been given in 

The fields ^apix) are expected to be related to the disorder fields fJ-aix) by the linear trans- 
formation 

/i«(x) = ^Cr/^p,(x), (33) 
fj 

where ClT £ C are coefficientJ§ to be investigated below, and p-independence is a consequence 
of permutational symmetry. The field fia carries the same representation of permutational 
symmetry as cTq (in particular, Yla=i t^ce = 0) but, as fiai3, it is non-local with respect to cTq. 
In other words, (Tq and fia are identical (dual) fields living in mutually non-local sectors of the 



^No confusion should arise with the OPE coefficients of ((30}. 
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3 3=(-+) 4=(++) 




Figure 2: The vacua of the Potts field theory for q = 2,3,4 are labeled by a = 1,. . . ,q and 
denoted by dots (for q = 4 we use the Ashkin- Teller notation cx = (01,02))- Fixing a specific 
vacuum (1 in this case) amounts to choose a representative kink field within each class fli (see 
the text). 



theory; in particular they share the same scaling dimension and the same OPE. Mutual non- 
locality reflects in the fact that, while (cJq) 7^ for J > Jc, 7^ for J < Jc- More precisely, 
in view of the coinciding scaling dimension, it is sufficient to adopt the same normalization of 
the fields to ensure that (ctq),/ = {na)j*, where J* is the dua0 of J. The duality extends to 
multi-point functions, in such a way that the spin correlators (jl2p can also be written as 

J*>Jc ■ (34) 

The rest of this section is devoted to investigate the relation ()33p and to determine the structure 
of the OPE fia • fJ-iB (or, equivalently, aa • (yp)- 

The q degenerate vacua of the Potts model above Jc can be associated to the vertices of 
an hypertetrahedron in g — 1 dimensions whose q{q — 1) oriented sides are associated to the 
kink fields ^ap- Permutational symmetry of the vacua allows to group these fields into classes 
fli, i = — 1, each containing q kink fields starting from different vacua, in such a way 

that choosing a vacuum amounts to select q — 1 kink fields, one from each class, starting from 
that vacuum and arriving at the other vacua (Fig. [2]). Ignoring structure constants, the OPE 
(31) has the form of a multiplication between elements of a finite group. Independence from 
the choice of the starting vacuum of the kinks ensures that the elements of this finite group 
are the classes /ij, i = l,...,g — 1, together with the topologically neutral class /. We denote 
then by Kg their fusion table as prescribed by (31), as well as the finite group of order q it 
defines. The symmetry also ensures that all the rows of the matrix Cq, can be obtained from 
the first by regular permutation^^. The relation (|33p (which we could equivalently write as 
fia = Ylp^c^"^ t^ptr) is effectively a sum over the q — 1 classes fli, as we now illustrate separately 
discussing the cases q = 2,3, 4. 

^''in the scaling limit we consider, J and J* are the points where the elementary excitations of the symmetric 
phase and those of the spontaneously broken phase have the same mass m; m = at the self-dual point Jc. 

^^Permutations which do not leave any element invariant. It is not difficult to realize that such permutations 
are elements of Kq. 
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Figure 3: Fusion tables Kq at g 
respectively. 



2,3,4. They correspond to the groups Z2, Z3 and 252) 



q=2. Permutational symmetry S'2 = Z2 implies the equivalence of the two kink fields ^12 and 
/i2i, which we collect in the class /ii; from (j32p we derive -R'2 = Z2 (see Fig. [3|). We can choose 



l^a — ^2 fXi, 







Wo 







1,2, 



with bJn 



Q2m/q^ in order to fulfill fia = and consistently derive 



(35) 



(36) 
(37) 



q=3. The OPE of the kink fields is equivalent to the fusion table of the classes jli,fl2 and /. 
Being Z3 the only discrete group of order three, full consistency requires = Z3, together with 
the identifications jli = {1M2, l-i23, M31}, 1^2 = {/^iSj ^21, /^32}- Notice that Aut(Z3) = Z2, and the 
non-trivial automorphisnj^i corresponds to the charge conjugation operator C, with Cfii = jl2- 



Usin 



35-. 



{x),a 



1, the charge conjugated operators realizing the Z3 OPE are identified 



witl^j = e^'^**'-^^/^ and fl2 = ^ 2-iTts{x)/3^ where s{x) is the dual color variable; using also 



IE 



13=1 



e 3 



■(s-a) 



one obtains 












^3 


^3 
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IJ,a = uj^^'fli +a;f/x2, 


Ca — 




^3 





W3-" 










^3 


■^3 





/ 



a = 1, 2, 3 , 



and then 



the relation a;? + uj^ + uj. 



-{"+/?) 



21 + Cf,Ha + 

-i - C^{lla + ^z?) + . . . 

= 0, a 7^ /3, is used. 



(38) 



(39) 
(40) 



q=4. The four-state Potts model can be seen as the case J = J4 of the Ashkin- Teller model 
defined by the Hamiltonian 



Hat 



{x,y) 



{J[Tl{x)Tl{y) +T2{x)T2{y)] + J4Tl{x)Ti{y)T2{x)T2{y)}, 



(41) 



^^Given a group G, : G — !■ G is an automorphism if (j){ab) = (j}{a)(j>{b),Wa,b G G. The set of all the 
automorphisms with natural composition as a product form a group called Aut(G). 
^^The basis jli, jl2 is that used in [TS] , 



11 



where Tj = ±1, i = 1,2, are Ising variables. Defining s = {ti,T2), cx = (ai,a2), with Oi = ±1, 
and 6s^a = ^Ti,ai^T2,a2i the Potts spin (fTU|) can be written as 

<7a = 4:6s,a. - 1 = aiTl + a2T2 + aiQ2Tir2. (42) 

The kink fields /Zq,^ interpolate between the four degenerate vacua of the two coupled Ising 
models (see e.g. [19]); the classes fli = {/xi2, /i2i, A*34, A'43} and /i2 = {mu, /i4i, /^23, A*32} are 
constructed in analogy to the case q = 2, the fields in jl-^ = {^13, /^ai, ^24, /^42} are instead 
obtained taking the OPE according to (31) (see also Fig. 2). We derive = 152 = ^2 x ^2, 
see Fig. [3l notice that Aut(D2) = S^. We can also take 



/ 

= aifli + a2fi2 + aia2/i3 , C, 








a\a2 


a2 


ai 





a2 


Qia2 


aia2 


a2 





ai 


02 


aia2 


ai 







a 







(43) 



from which we obtain 

^loc ■ ^i» = 3I + 2C^Hc. + ■■■ , (44) 

(45) 

It is interesting to remark some formal properties emerging from this analysis. We see that, 
by construction, Ka at q = 2, 3, 4 is a finite abelian group of order o, i.e. by Cayley theorem a 
regular abelia. subpouB of K, m.s. also be ™t under permutations of the , - 1 
classes /2j, an operation which corresponds to fix one vacuum and permute the remaining q — 1- 
Formally this amounts to write Aut{Kq) = and we expect the full symmetry group of the 

theory to be realized aJ^ 

Sg = Kg>4 . (46) 

This in turn implies the possibility of writing Sq as a semidirect product of abelian subgroups 
of the form 

5g = Kg X Kq.i X • • • X ^2, K2 = Z2, (47) 

a property which is equivalent to the solvability of the permutational group. More preciselj0, 
the solvability of Sq would imply the existence of the factorization ()47p . as indeed remarkably 
happens at q = 2,3,4, with K2 = ^2i K3 = Z3, K4 = D2 and Aut(Z3) = Z2, Aut(L'2) = S3. 
It is well known |21] . however, that for q > 4 Sq possesses no abelian normal subgroup, making 
impossible in particular to realize the condition (|46p . Then Sq is not solvable for any q > 4, a 
circumstance which is interesting to compare with the fact that = 4 is also the upper bound 



^''The classes p,i are associated to the regular permutations ivi, i = 1, . . . , of Sq as jii = {/ii^.(i), ^^^^(g)}. 
Without loss of generality one can assume 7ri(l)=i + l. 



"'The presence of the semidirect product xi is due to the fact that Sq-i is not a normal subgroup of S, 
^We thank C 
be found in [20| 



^®We thank C. Casolo for this observation. Interesting remarks about solvable groups and lattice duality can 
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for the existence of the Potts field theory in the two-dimensional case, i.e. the only case for 
which kink fields exist and the above construction is possible. 

A remarkable feature appearing from (f37|) . ([I0|) and psj) is that fia • fJ-is^a is identical at 
q = 2,3,4 (recall that /^i + = at (7 = 2). It is then absolutely natural to assume that this 
form actually holds unchanged for continuous values of q, and to write 

(^a- (^13 = -I - Ci,{aa+(7is) + ... , Oy^/S, (48) 

where the dots correspond to less relevant fields and we switched to the equivalent expression 
in terms of the spin fields for a reason to be made immediately clear. On the other hand, the 
complementary relation 

(Ja ■ (7a = qi I + q2Cf, (Tq H (49) 

follows observing that ([IT]) and give 

= CTq • ^<T/3 = <Ta • (To + ^ (Ta- O-fS 

= CTa • (Ta + ^ [-/ - C^{aa + (7 p) + . . .] ] 

(fl9l) is of course consistent with (|36]l . (f39l) and (jM]). While the disorder fields fJ-a{x) are specific 
of the two-dimensional case, the spin fields aa{x) are well defined in any dimension. It is then 
quite obvious to expect that (j^Hj) and (09]) hold for real values oi q < qc in any dimension. 

The linear relation (jlip among the spin fields induces a relation less direct than usual between 
the OPE coefficients and the structure constants appearing in the three-point functions. In 
general, the structure constants Cijk are defined by the critical correlators 

C ■ 

{Mxi)Aj{x2)Ak{x3)) = x.+t,-x, x,+x,-x, ; (50) 

■^12 '''13 "^23 



taking the limits xi2 — )• and 2:23 — )• and using (j30p . Cijk is expressed in term of the OPE 
coefficients as 



Xm—Xi.^ 



We then find 



Ca^a^a^ = q2qiC^, (52) 

Ca-aaaap = Ca-^o-is<Ja — Ccrpcracra ~ ~q2C^, (53) 
Caaa,sa^ = 2C^ , (54) 



with different indices denoting different colors. 
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4 Duality relations 



Equations (|33p and (j34p imply a linear relation between the spin correlators (|12|) in the symmetric 
phase and the correlators (j27p of kink fields in the broken phase. This duality takes the form 

/3l,...,/9n 

nc(/3)-l 

= E' ( n Dt-:±Gp,...p^{xr,...,Xn), (55) 

/3i,..,/3„ i=l 

where the primed sum runs over all choices of (3 = {/3i, . . . , /?„} which are inequivalent under 
permutations, ndP) is the number of different colors in /3, and the factors have been extracted 
for later convenience. The task is that of determining the coefficients Da\'.'.'.aZ: for continuous 
values of q; of course it is sufficient to consider a set of linearly independent spin correlators. 
We will discuss explicitly this problem up to the first case with F„ > 1, i.e. n = 4. 
n=l. The trivial identity (fJo) = (fia/^) = simply reflects the fact that we consider spin 
correlators in the symmetric phasfl , and that n^p is a kink field. 
n=2. F2 = 1 and we consider 

{(7a{xi)aaix2)) = QlD'^^ (/^q/? (a^l )/U/3q (^2)) . (56) 

We can take the limit X12 — )■ on both sides and equate the coefficients of the leading singularity 
Xi2^^'^ , which are obtained using (P9|) and (j32]) . This immediately yields Daa = 1- 
n=3. F3 = 1 and we consider 

{craixi)aa{x2)(7a{x3)) = qiq2D^al (^^1 )M/37 (^^2)^70 (2:3)) • (57) 

We can again use the OPE's to take the limit xi2 — t- on both sides and reducj^ to (|55p : this 
leads to DaaZ = 1- 

n=4. We will consider the F4^ = 4 linearly independent spin correlators ()19H22p . which expand 
as 

((7ai(xi)ac,2(x2)cTQ3(x3)fJ„4(x4)) = QmqsD^i^LaaAf^'^pi^l) f^l3^i^2) fJ-^si^s) fiSa{x4:)) 

+9l92-D"fjf<^3,,^ (Ma/3(a;i)/i/37(2^2)/i7/3(x3)/i/3a(x4)) 

{f^ap{xi)fli3c,{x2)f^a-,(.X3)fl^a{x4)) (58) 
{tJ'al3{xi)tJ-l3a{x2)fJ.al3{x3)fIi3aixi)) ■ 



We start again equating the coefficients of the short distance singularities on the two sides. 
However, an important difference with the cases n < 4 now appears. Indeed, while the OPE 

^^Since this is understood, here and in the following we omit the subscripts J < Jc for spin correlators and 
J > for kink field correlators. 

^^Notice that the OPE on the l.h.s. apparently produces singularities from the S'q-invariant operators Ok in / 
which do not arise in the r.h.s., where 0^7. Everything is consistent, however, since {Okcra) = by symmetry. 
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Figure 4: Pictorial representation of the four-point correlation functions of kink fields 
^^i/32fe/34(^i' • • • ' ^4)- The two pinchings xi — )• X2 and X2 — )• X3 used for the study of the 
duality relation (j58|) are indicated. 



in the l.h.s. can be taken for any pair of points Xj and xj, the kink nature of the fields on the 
r.h.s. allows us to use (1320 only for adjacent fields. In other words, we can compare only the 
singularities arising for X12 — )■ and X23 — )• (see Fig. HI 2:34 — )• and X14 — >• add nothing 
new). Equating the coefficients of the singularities in the "neutral" and "charged" channels, and 
using the duality relations already obtained for n = 2, 3, leads to the following sets of equations 



^aaaa ^aaaa 

09 = 03 L*"^'^'' + 



1 - gmJ^afiafi + 9lJJal3al3 



(59) 



1 = 

2 = 
k 91 



(60) 



where different indices denote different colors. Notice that the method produces four equations 
for each correlation function (two per pinching and per channel), but only three turn out to be 
independent. Since each duality relation (j58|) involves four coefficients, the above equations are 
not sufficient to fix everything. We now show how the duality for the correlators Gaaaa, Gaapp-, 
Gap Pa can be completely determined exploiting also the relations (p3]) . (jMj) . 

The matrices Ca defined in (|33p are hermitian, due to the presence of the antilinear charge 
conjugation operator C, with = /, C^p^C = fi^p, CfiaC = fJ-a and CCa^C = {C^)* . They also 
satisfy Ca'' = 0, p = 1, . . . , and Ca = 0. Other properties follow requiring the consistency 
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of the OPE's for the kink fields and for the dual spin /ia. For example, 

(61) 

cr,aj 

E cr) + E ( E ^r^^r) Kf^p^ + ...]; (62) 

comparing with (j49p we obtain the matrix relation 

Cl = qiI + q2Co.. (63) 

Analogously, from ()48p we derive 

C„C^ = -/ - (C„ + C;3) , a7^/3, (64) 

which implies in particular [Cq,,C^] = 0. Hence, the set of q hermitian matrices can be 
simultaneously diagonalized by a unitary transformation U and put in the form 

C^ridiag = (9'^«p-l)V, (65) 

which follows from the observation that the Cq's are traceless, sum to zero and, due to (|63p . 
have —1 and qi as only eigenvalues. It is also simple to check that 

(/^a/^o) = ~^^^a = 9l (/"a/3/"/3Q:) , (66) 

1 „ 

(/Xa/Xa/ia) = -TrC^ (/^a/3/^/37/"7a) = 9291 (/"a/3 /U/37/^7«) , (67) 

in agreement with (j56p . ()57p . Starting from the diagonal form (j65p and the existence of the 
unitary matrix U it is possible to show (see Appendix C) that 

CP," = e*(^"--^-'') - 5p„, (68) 

where the phases ipap must satisfy the equatio: 

1 ^ 

= § (69) 

We can now notice that (f33]) . (p^ and ([55]) imply 

1l^a^a2a3a4, = E ^"i ^02^03 ^'af ; (70) 

using IC^""! = I - 5pa and hermiticity Cq^ = {C^")* we obtain L>Sa2a = D'^^jfp = D'^p'fa = 1' 
independently of the phases i^ap- With this information, ()59p and ()60p determine the remaining 



^^The properties of matrices Ca we obtain in this section do not refer to any specific value of q. Of course they 
are satisfied by the matrices that we aheady determined in the previous section for q = 2, 3, 4. 



16 




Figure 5: Four-point boundary correlations on a simply connected domain L. Left: Clusters 
connecting xi to X3 and X2 to necessarily cross and cannot be distinct. Right: Dashed 
lines represent allowed clusters in the FK representation on the dual lattice £*; continuous lines 
correspond to clusters on C connecting all the four points. In the continuum limit C, C* and L 
coincide. 



coefficients, giving 

Gaaaa = 9l92'73 Ga(3yS + QlQ2 {Oap^p + GaPa-y^ + Ql Gapap, (71) 

GaaPP = — 9293 Ga/i-yS — 92 Ga/i-y/B + QlQ2 Gaf^a-y + Ql Ga/3al3, (72) 
Gappa = —<i2<i3 GafS-yS + <il<i2 GalS-y/B — Q2 Gaf^a-y + Ql Ga/3al3, (73) 

with different indices denoting different colors. 

The problem with the remaining correlator Gapap is that (I70p does not help, because the 
phases do not cancel, and we are left with the three equations coming from the OPE for four 
unknowns. 



5 The boundary case 

Since the OPE's for the kink and spin fields reflect local properties of the field theory, the 
duality relations obtained in the previous section hold true also in the case in which the points 
xi, . . . ,Xn in (|55|) . instead of being located on the infinite plane, lie inside a simply connected 
domain L C M^. Actually, the duality relations continue to hold also in the case the points 
located on the boundary of L, simply because the OPE's ()32p . ()48p and ()49p . 
whose structure is completely determined by the symmetry, can be used also for points on the 
boundary, provided bulk OPE coefficients and scaling dimensions are replaced by boundary OPE 
coefficients and scaling dimensions. It is not difficult to see, however, that having the points on 
the boundary rather than in the bulk may reduce the number of linearly independent correlation 
functions. The first interesting case, that we now discuss, arises for n = 4. 

Let us order the points xi, . . . , X4 on the boundary as shown in Fig. [SJ Since in this boundary 
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case a cluster containing xi and X3 must necessarily cross a cluster containing X2 and X4, the 
probability Pababi^i, x^, X4) that these two pairs of points belong to two different clusters 
necessarily vanishes. This topological constraint reduces to three the number of linearly inde- 
pendent boundary spin correlators (we denote them with a superscript B); indeed, inverting 



(1T9|)-(|22]) and setting Pabab = give^ 

Qlig"^ - 3g + l)G^i3al3 - (2? - ^)Gaaaa + Ql{Gaal3l3 + Gai3l3a) = • (74) 

Using the duality relations ([7T]) . (j72|) . ([73]) and the OPE equations (j59|) for Ga^ap, can be 
rewritten as 

[1 + gi {q' -3q+ 1)02^] {G%^^ + G%^, - G^, " G^^.p) = , (75) 
or, in view of the reduction in the number of independent correlators , as the linear relatior 

^a/3a/3 + '5f/3^5 = Gf^Q,^ + Gf^^^. (76) 

Using this equation to eliminate G^^^p in ([7T|). ([72]) and ([73]) . the duality for boundary corre- 
lators is fully determined as 

Gaaaa = <?l(9293 " 1)Gq^^5 + QliGap-ylS + ^f/jQ,^), (77) 
Gaal3l3 = (292 — Ql)Gal3'y5 + GajS-^jS + QlGapa'yJ (78) 
Cf^/Jo = (292 - Ql)Gal3^s + QlGaf3^i3 + G^l3a^- (79) 

Substituting in p^ - (|22]) with Pafeafo = one also obtains for the boundary connectivities the 
simple relations 

(80) 
(81) 
(82) 

Equation (j80p can be interpreted as follows in the language of lattice duality. If we interpret 
the insertion of a kink field on the boundary as creating a domain wall along the boundary of 
the dual lattice C*, the correlator G*^^^^ corresponds to a partition of the boundary into four 
regions with different colors and will receive contributions only from graphs on C* without FK 
clusters connecting different regions. Equation (j80p then means that these graphs are in one-to- 
one correspondence with the graphs on C in which the four boundary points all belong to the 
same FK cluster (see Fig. [5])). Similar reasonings can be used for (181 p and ()82p . 



pB 
^ aaaa 


- 


1 


pB 
^ aabb 


_ fiB 


f^B 


pB 
^ abba 


r<B 


r<B 



■^"DifFerent greek indices in eqs. H74H82p denote different colors. 

■^^ Duality for boundary correlators of the g-state Potts model on the lattice was studied in [221 1231 [24] . where 
dual partition functions with domain wall boundary conditions correspond to our kink field boundary correlators 
Gaia2,...\ with this identification, the relation (|76|l is contained in [21]. Potts partition functions on a non-simply 
connected domain have been studied in [25]. An early investigation of Potts correlation functions is in [26] . 
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a 



t= 



y 



r=-oo 



Figure 6: The elastic kink- kink scattering amplitudes Sa^{s) (left) are related by the LSZ 
formalism to the four-point correlation functions G^/^-yS of the asymptotic fields /x^^(x) (right). 



We conclude this section observing that (|76p leads to a linear relation among the elastic kink- 
kink scattering amplitudes in the (1 -|- l)-dimensional field theory associated by Wick rotation to 
the Euclidean theory on the plane. In this case, indeed, the boundary fields /i^^(ic) entering ([76]) 
can be interpreted as the asymptotic fields which create a kink excitation Ka/s at time t — )• ±00 
(Fig. E]). The kink- kink elastic amplitudes 

^^(S) = ont{Ka5iPA)Ksy{p3)\K^M)Kp^{p2))in , (83) 

where the square of the center of mass energy s = {pi +P2)fi is the only relativistic invariant for 
this (1 -|- l)-dimensional process, can be written within the LSZ formalism (see e.g. [27j) as 

5f^(s) = lim lim f dxi... ( dx4 e^^'^-^^gips-^Sg-ipi-xig-ipa-xa 

il,i2— oo 43,44— >+oo J J 

^ ^ ^ ^ (^f^(xi,ii)/i^^(x2,t2)/^:f5(a;3,t3)Aifa(2;4,t4)> , 



(84) 



where the integrals are taken along the one-dimensional space coordinate and dt is defined by 
A^tB = A{dtB) - {dtA)B. Equation ([TeD then leads to 

S^'^{s) + S^iis) = S^/^is) + S^l{s) , (85) 

where different indices denote different colors. This relation was already found in [5] as a 
byproduct of the Yang-Baxter equations, i.e. of the integrability of the scaling limit of the q- 
state Potts model. Here integrability does not appear to play a role, and then (|85p should hold 
also for non- minimal, non-integrable realizations of the symmetry, if they exist. Ultimately, the 
vanishing of the probability that the trajectories of two particles do not cross in the (1 -|- 1)- 
dimensional space-time relates amplitudes which would be independent on the basis of color 
symmetry alone. 



19 



6 Conclusion 



In this paper we considered the g-color Potts field theory characterized by invariance under 
color permutations and having the spin fields daix) as fundamental fields. We discussed how 
the theory can be given a meaning for real positive values q < Qc for the purpose of describing 
the scaling limit associated to the second order phase transition of the random cluster model, 
consistently with the expectations coming from the Fortuin-Kasteleyn correspondence. In par- 
ticular we showed that, although the number of spin fields is badly defined for q non-integer, the 
number of linearly independent n-point spin correlation functions is a g-independent integer F„ 
coinciding with the number of partitions of n elements into subsets each containing more than 
one element. These n-point spin correlators, which in turn determine the connectivities in the 
random cluster model, contain ric < n different color indices. Hence, they enjoy an ordinary 
definition for integer values oi q > Uc, but formal use of the symmetry unambiguously prescribes 
their analytic continuation to real values of q. The mechanism allowing the analytic continuation 
is exemplified by the form of OPE that we obtained for the spin fields (equations (j48|) . (|49|) ): 
given two colors, the symmetry only discriminates whether they are equal or different, and q 
enters as a parameter which can take real values. 

We also discussed how the purely magnetic properties of the Potts model, which are well 
defined for integer values of q, are described by the sector of the theory in which the only 
spin variables aa playing a role are those with a = 1, . . . ,g. The number Mn{q) of linearly 
independent n-point correlators of such variables, which is smaller than F„ for n large enough, 
has been determined for q = 3, 4. 

For the two-dimensional case, we studied the duality relations between the spin correlators 
in the symmetric phase and the kink field correlators in the spontaneously broken phase, both 
for bulk and boundary correlators, exploiting the OPE as only input. For the case of four-point 
correlators, the first with F„ > 1, our method gives the duality relations for three of the four 
independent bulk correlators, while it completely fixes the duality for independent boundary 
correlators, whose number is reduced to three for topological reasons. The ultimate reason for 
our inability to fix completely the bulk duality is that the OPE ()32p for kink fields determines 
the fusion of adjacent fields only: while this is the only possibility on the boundary, in the bulk 
also the fusion of non-adjacent fields (i.e. xi — )• X3 in Fig. H]) is allowed. Clearly this point 
requires further investigation. 

Finally, we showed how the constraints that topology imposes on boundary correlators leads, 
through the LSZ formalism, to relations among kink scattering amplitudes in (l + l)-dimensional 
space-time. These relations, already observed within a framework based on integrability, appear 
in our derivation as a generic feature of two-dimensional field theories with spontaneously broken 
Sq symmetry. 

Acknowledgments. J.V. thanks A. De Luca and F. Mancarella for discussions. Work sup- 
ported in part by the MIUR Grant 2007JHLPEZ. 
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Appendix A 



Given a set S of n elements the number of partitions of the elements of S is the Bell number 
Bn (see Table [T] for n < 10). The most straightforward way to compute the Bell numbers is 
through the recursive relation 



n-l , _ 
fc=0 



which is easily proved observing that we can fix one element of the set S and consider the 
partitions in which this element appears with k other elements, A; = 0, . . . , n — 1. The number 
of such partitions will be 

^ ^ Bn-i-k ; (87) 

then summing over k and using {^^^) = (^"7-fe) obtain (f86]l . 
We similarly define the numbers F„ as 



^ \k 

k=0 ^ ' 

Fn is the number of partitions of S whose blocks contain at least two elements. The proof is 
again elementary (see e.g. |28]). We divide all the Bn partitions of S into those containing 
exactly k = 0,1, . . . ,n isolated points, and then we take partitions of the remaining n — k points 
in such way that no point is isolated. It is clear that we end up with (f88]l . 

Recalling the combinatorial identity ("^^) = (^) + (;,"J, expression (|88p implies 

Br,+i-Bn = ^(''^Fk+i. (89) 
k=o ^ ^ 

On the other hand, (f86]l gives i?n+i = Sfc=o (fc)-^*:' ^^'^ using (f88|) we derive 

n / 

X - / n 



fc=0 
221 



^ (Fk+i + Fk - Bk) = Vn>0. (90) 



By induction we finally obtai: 

Bn = Fn + Fn+l . (91) 

The number of /c-partitions (partitions into k non-empty subsets) of a set of n elements is 
the Stirling number S{n,k). It satisfies Bn = Z^fc=i 'S'(n., A;) from its definition, as well as the 
recursive equation 

S{n,k) = kS{n-l,k) + S{n-l,k-l), n>k, k>l, (92) 



^■^Alternatively, one can recover (|9ip introducing the exponential generating functions B{x) = Y2„ ~ ®° ~^ 

and J-{x) = 7^2;" = e'' the result then follows from B{x) = J-{x) + T'{x). We thank the referee for 

this observation. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Fn 





1 


1 


4 


11 


41 


162 


715 


3425 


17722 


S{n,l) 





1 


1 


1 


1 


1 


1 


1 


1 


1 


S{n,2) 











3 


10 


25 


56 


119 


246 


501 


5(n,3) 

















15 


105 


490 


1918 


6825 


S{n,4) 























105 


1260 


9450 


S{7i, 5) 





























945 



Table 2: The Fn independent n-point spin correlators without isolated indices decompose into 
subsets containing S{n,k) correlators with indices of k different colors. 

from the fact that we can obtain a A:-partition of the set {xi, . . . , Xn} adding x„ to one of the k 
blocks of a /c-partition of the elements {xi, . . . , Xn-i}, or joining Xn as a single block to a (A; — 1)- 
partition of {xi, . . . ,Xn-i} ■ The exponential generating function Sk{x) = l^n>fc 5'(n, /c)^ 
satisfies = kSk{x) + Sk-i{x), and is given by Sk{x) = ^(e^ — 1)'^- We also introduce the 

generalized Stirling number S{n, k) as the number of fc-partitions of a set of n elements whose 
blocks contain at least two elements (non-singleton A:-partition) ; the relation 

n-l 

Fn = Y,~S{n,k) (93) 
fc=i 

then expresses the decomposition of the total number of independent n-point spin correlation 
functions (which we take without isolated indices) into subsets with indices of k different colors. 
Non-singleton ^-partitions of the set {xi, . . . , x„} are obtained adding Xn to one of the k blocks of 
a non-singleton fc-partition of {xi, . . . , x^-i}, or by joining the block {x^, x^}, for j = 1, . . . , n— 1, 
to a non-singleton {k — l)-partition of {xi, . . . , Xj_i, x^+i, . . . , x^-i}. We have then 

5(n, k) = kS{n -l,k) + {n- l)S{n - 2, - 1), n>k,k>l. (94) 

The exponential generating function is 5fc(x) = -^(e^ — 1 — x)'^' and solves 5^(x) = kSk{x) + 
xSk-i{x). The first few S{n,k) are collected in Table [2j 

Appendix B 

In this appendix we determine the number Mn{q) of S'g-inequivalent, linearly independent n- 
point spin correlators (jl2p which determine magnetic correlations in the Potts model at g = 
2,3,4. We exploit the fact, discussed in Section 3, that at q = 2,3,4 the symmetric group 
factorizes as 

= Kg X , (95) 

with K/i = D2, = Z3, K2 = ^2, and that the Potts model is described by g — 1 indepen- 
dent spin variables ti, . . . charged under the abelian group Kq. Non- vanishing correlation 
functions are neutral under Kq and invariant under permutations of the q — 1 operators ij. 
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At g = 2, the only independent variable is ti with charge +1 under Z2. The neutrality 
condition for the n-point correlation functions is 



n = mod 2, (96) 

giving M2fc(2) = 1, M2fc+i(2) = 0. 

At g = 3, the independent spin variables are ti and t2 with Z3 charge +1 and —1, respectively. 
Given a ra-point correlation function containing ni variables ti and n2 variables t2 we require 

ni + n2 = n, (97) 
ni — = mod 3, (98) 

or equivalently ni+n = mod 3, with ni = 0, . . . , n. Assigned the couple of integers {^i, 722} = 
{ni,n — rii} satisfying the constraint ()98p . the total number of distinct correlation functions 
we can construct is the binomial coefficient (^^) . The charge conjugation operation, however, 
exchanges rii with 722, and correlation functions obtained by ni — )• n — ni are equal; notice 
indeed that (") = ("„). It followfl 

\ni/ \n—n\/ ^ — ' 

t {:)■ (») 

ni+n=0 
mod 3 

The elements of the sequence (1990 (see the first few of them in Table [I]) coincide with the 
Jacobsthal numbers [29] and satisfy the recursive relation 

M„+i(3) = Af„(3) + 2M„_i(3) , (100) 

with Ml (3) = 0, M2(3) = 1. We will not prove (jlOOp directly but we will justify it through 
the following observation. Consider an hypertetrahedron with q vertices labeled by the numbers 
!,...,(?. The number y^") of closed n-step paths starting from a given vertex, say 1, of the 
hypertetrahedron satisfies the recursive relation (see Appendix C) 

y(") = (g-2)y("-i) + (g-l)y(-2), (101) 

with y^^^ = 0, 2/^^^ = g — 1. The closed n-step paths 

in (jlOip are considered distinct even 
when they differ by a permutation vr G Sq-i of the q — 1 vertices 2, . . . , g. At g = 3 there is 
only one possible permutation vr, and it exchanges the vertices 2 and 3. The application of vr 
to a path 7^") generates the path reflected along the symmetry axis containing the vertex 1 of 
an equilateral triangle. The number of closed paths inequivalent under permutations at g = 3 
is then just half of the total number of closed paths y*-"-*, and in particular satisfies (jlOip with 
(7 = 3. Finally notice that closed n-step paths 7^"^ inequivalent under permutations are in one 



^''if n is even and ni = n/2 the factor 1/2 in (|99|) avoids the double counting of the correlation functions 
obtained exchanging in block the positions of the ni operators t\ with the 712 = n\ operators t2. 
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to one correspondence with independent n-point kink fields correlation functiona^i (j27p . whose 
number is M„(g). Consistency requires that M„(3) satisfies the recursive equation (|1U1|) with 
q = 3, which indeed coincides with (jlOOp . 

For q = 4 we consider correlation functions of the three variables ti , t2 and t-s with D2 = 
7L2 X Z2 charges (1, 0), (0, 1) and (1, 1), respectively. A non-vanishing n-point correlation function 
with m variables ii, n2 variables t2 and ns variables tz satisfies 

n\^n2^n-i = n, (102) 
{ni + n-s, n2 + n-s) = (0,0) mod 2 , (103) 

or, more symmetrically, rii = n mod 2, i = 1,2,3. The number of distinct n-point correlation 
functions associated to the solution {ni,n2,n3} of (|102p and (jl03p is "jingi ■ Correlation 
functions obtained by permutations of the integers nj are identified and we must therefore 
choose a definite order for them, for example ni < n2 < n^; nj = 0, . . . ,n. When two positive 
integers ni and n2 coincide, the two correlation functions obtained by exchanging in block the 
positions of the rii operators ti with the n2 operators t2 are also equal and counted twice among 
the ^^i^ji^gi correlation functions. Similarly, if ni = n2 = permutational symmetry does not 
distinguish among the 3! correlation functions obtained exchanging in block the positions of the 
Hi variables tj for i = 1, 2, 3. The final result is then 

„,+irn3=n"i'"2!n3! ne(ni, n2, n3)! 

ni=n mod 2 

where the nj are ordered, n, < nj+i, and ne(ni, n2, ns) is the number of non-zero equal integers 
in the tern {ni, 712, ""-a}- The integer sequence (jl04p (see Tabled]) is also known [30], and is 
solution of the recursive equation 

M„+i(4) = 2M„(4) + 3M„_i(4) - 1 , (105) 

with Ml (4) = 0, M2(4) = 1. Again we will not prove (|105p directly, but will explain why 
a proper counting of closed paths ^^'^'^ inequivalent under permutations of the vertices on a 
tetrahedron is obtained subtracting 1 to the r.h.s. of (|10ip with g = 4. It is convenient to 
represent 7*^*^) as the sequence 7^"^ = {1,^2, . . . iVn, 1} with Vi = 1, ... ,4 and Vi 7^ Uj+i, which 
is also a particular coloration of n points on a circle. The way an (n + l)-step closed path 
^('i+i) = {l,f2, . . . ,Vn,Vn+i, 1} is Constructed by adding a new point Vn+i and taking care of 
permutational symmetry is the followinj^. First we identify the largest f j, i = 1, . . . , n, and then 
generate all the paths 7^"+^) for which Vn+i / 1, with Vn+i = 2, . . . ,min{max{uj} + 1,4}. 



^^Any path 7*-"^ can be represented as the sequence of n + 1 vertices {1,V2, ■ ■ . ,Vn,l} with Vi 7^ u^+i and 
Vi = 1 . . .q. The associated kink fields correlation function is {fi\v2{xi) ■ ■ ■ fJ-vnii^n)) ■ Alternatively 7'"-* can be 
thought as a particular coloration of n points on a circle. See again Appendix C. 

■^^Notice that there are two cases corresponding to i;„ = 1 or w„ 7^ 1. In the first case the new point Vn+i is 
added to some closed (n — l)-step path 7("-i'; in the second case the new point is added to an (n — l)-step open 
path. The number of (n — l)-step open paths is however equal to the number of closed n-step paths 7'"' 
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The recursion (jlOip fails only if max{uj} = 2 and then Vn = 1,2. In this case the closed paths 
^("+1) with Vn+i = 3, 4 are identified by permutational symmetry and the choice fn+i = 4 must 
be discarded; this leads to (jlOSp . 

We conclude noticing that n distinct vertices are sufficient (and necessary) to enumerate 
all closed n-step paths inequivalent under permutations of the vertices on an hypertetrahedron. 
In particular the number Mn{q) of closed n-step paths inequivalent under permutations on an 
hypertetrahedron with q vertices is constant for any q > n, and we have already shown in Section 
2.3 that, when no restriction is assumed on the number of available vertices so that all closed 
inequivalent paths are counted, this number coincides with Fn'. 

Mn{q)=Fn for q>n. (106) 



Appendix C 

Given the diagonal form (j65p and the unitary matrix [7, we have 



^ pa'-^ ctcF '-'pa- 



Requiring Ca'' = for p, a = 1 . . . g, gives 



ap 



V9 



(107) 
(108) 

(109) 



the equation (|69p for the phases is then the self-consistent condition of unitarity of the matrix 
U. Substituting (fT09]l back into (fT08]) we obtain The matrices C„ in ([35]) and ([MI) for 

q = 2,3 are reproduced by the solution 



2tt 

^ap = ± — ap 



(110) 

of the phase equation ([69]) . For q = 4, the matrix (jl3]) corresponds instead^^l to the solution 

ipocp = ±7r(ai/)i + a2P2), (111) 

with a = (01,02), p = {pi,P2), ai,pi = 1,2. 

We conclude this appendix giving a simple geometrical interpretation of the relation (j63p . 
Without loss of generality we can choose the phases so that one of the matrices Cq, say Cq, is 
real, i.e. 

/O 1 ••• 1\ 
1 ••• 1 



V 1 



(112) 



0/ 



^''Of course (|110|) solves (|69|) also for g = 4. This solution, however, would lead to a matrix Ca associated to 
K4, = Z4, inconsistent with our general discussion of section 3 (in particulax, Aut(Z4) = Z2). 
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This is the adjacency matrix of a fully connected graph with q sites, i.e. the projection on the 
plane of a hypertetrahedron with q vertices. We can then use the formula (|63|) to compute 
powers of Cq and obtain 

= y (")/ + (113) 
with the integers x^") and y^"^^ satisfying the recursive equations 

= (114) 

=g2x(")+y("). (115) 

It is now simple to realize that y^"^ and x*-"^ are, respectively, the number of closed paths of 
length n starting from a vertex vi of the hypertetrahedron and the number of open paths of 
length n from vi to Vk- Indeed, if we consider such an open path and we add the link {vk,vi) 
we obtain a closed path of length n + 1; however, by permutational symmetry, all the locations 
of Vk 7^ vi are equivalent, and we obtain ()114p . Suppose instead to remove from the original 
open path the last link {vj,Vk)', this gives a path of length n — 1 from vi to vj that can be open 
with multiplicity q2 {j ^ 1,^), or closed with multiplicity one (j = 1), reproducing (jllSp . The 
recursions (jll4p and (jllSp give for y^") equation (jlOip . whose solution is 

yW = i[(_l)ng, + gn], (116) 

q 

which is q~^ times the chromatic polynomial^ nn,^ (q) for the the cyclic graph C„, i.e. the graph 
obtained putting n points on a circle. We conclude that the decompositions over kink fields 
of the correlator (/Xq,(xi);Uq,(x2) . . •/io(xn)) can be associated to closed n-step paths on a fully 
connected graph with q vertices, or to colorations of a ring of n points with q colors. 
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